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Abstract 

' In this work we establish some new estimates for layer potentials of the acoustic 

. wave equation in the time domain, and for their associated retarded integral op- 

fH I erators. These estimates are proven using time-domain estimates based on theory 

' of evolution equations and improve known estimates that use the Laplace transform. 
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G^ : 1 Introduction 
a^ 
m 

Tj- ■ In this paper we prove some new bounds for the (two and three dimensional) time domain 
Q . acoustic wave equation layer potentials and their related boundary integral operators. 

In 1986, Alain Bamberger and Tuong Ha-Duong published two articles (references [2] 
and [3j) on retarded integral equations for wave propagation. These seminal papers estab- 
lished much of what is known today about retarded layer potentials, proving continuity of 
layer potentials and their associated integral operators as well as invertibility properties of 
some relevant integral operators. The analysis of both papers has two key ingredients: (a) 
the time variable is dealt with by using a Laplace transform; (b) estimates in the Laplace 
domain are proved using variational techniques in free space, very much in the spirit of 
[T5] (see also [E]). Even if the results in |2] and [3] are given only for the three dimen- 
sional case (retarded operators with no memory), because of the way the analysis is given, 
all results can be easily generalized to any space dimension. An additional aspect that is 
relevant in ^ and [3] is the justification of time-and-space Galerkin discretization of some 
associated retarded boundary integral equations, a result that sparked intense activity in 
the French numerical analysis community on integral methods for acoustic, electromag- 
netic and elastic waves in the time domain. Not surprisingly, when Lubich's convolution 
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quadrature techniques started to be applied to retarded boundary integral equations (this 
happened in |I3]), the key results of Bamberger and Ha-Duong were instrumental in prov- 
ing convergence estimates for a method that relies heavily on the Laplace transform of 
the symbol of the operator, even though it is a marching-on-in-time scheme. The rele- 
vance of having precise bounds in the Laplace domain for numerical analysis purposes has 
also been expanded in more recent work at the abstract level (with the recent analysis of 
RK-CQ schemes in [5j and [6]) and with applications to the wave equation at different 
stages of discretization ([12], [8]) 

In this paper we advance in the project of developing the theory of retarded layer 
potentials with a view on creating a systematic approach to the analysis of CQ-BEM 
(Convolution Quadrature in time and Boundary Element Methods in space) for scattering 
problems. As opposed to most existing analytical approaches -while partially following 
the approach of [E]-, we will use purely time-domain techniques, inherently based on 
groups of isometries associated to unbounded operators and on how they can be used 
to treat initial value problems for differential equations of the second order in Hilbert 
spaces. We will show how to identify both surface layer potentials with solutions of wave 
equations with homogeneous initial conditions, homogeneous Dirichlet conditions on a 
distant boundary and non-homogeneous transmission conditions on the surface where the 
potentials are defined. This identification will hold true for a limited time-interval, and 
a different dynamic equation (with a new cut-off boundary placed farther away from the 
original surface) has to be dealt with for larger time intervals. In its turn, this will make 
us be very careful with dependence of constants in all bounds with respect to the (size of 
the) domain. Bounds for the solution of the associated evolution equations will depend 
on quite general results for non-homogeneous initial value problems. A delicate point 
will be proving that the strong solutions of these truncated (in time and space) problems 
coincides with the weak distributional definitions of the layer potentials. Since the type of 
results we will be using are not common knowledge for persons who might be interested 
in this work, and due to the fact that the kind of bounds we need are not standard in the 
theory of Co-semigroups (and, as such, cannot be located in the best known references on 
the subject), we will give a self-contained exposition of the theory as we need it, based 
on the simple idea of separation of variables, the Duhamel principle, and very careful 
handling of orthogonal-series-valued functions. 

From the point of view of what we obtain, let us emphasize that all bounds improve 
results that can be proved by estimates that use the Laplace transform. Improvement 
happens in reduced regularity requirements and in slower growth of constants as a function 
of time. This goes in addition to our overall aim of widening the toolbox for analysis of 
time-domain boundary integral equations, which we hope will be highly beneficial for 
analysis of novel discretization techniques for them. 

Although results will be stated and proved for the acoustic wave equation (in any 
dimension larger than one), all results hold verbatim for linear elastic waves, as can be 
easily seen from how the analysis uses a very limited set of tools that are valid for both 
families of wave propagation problems. Extension to Maxwell equations is likely to be, 
however, more involved. 

The paper is structured as follows. Retarded layer potentials and their associated 
integral operators are introduced in Section [21 first formally in their strong integral forms 
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and as solutions of transmission problems, and then rigorously through their Laplace 
transforms. Section [3] contains the statements of the two mains results of this paper, one 
concerning the single layer potential and the other concerning the double layer potential. 
Sections S] and \5\ contain the proofs of Theorems 13.11 and 13.21 respectively. In Section 
E] we use the same kind of techniques to produce two more results, much in the same 
spirit, concerning the exterior Steklov-Poincare (Dirichlet-to-Neumann and Neumann-to- 
Dirichlet) operators. In Section [7] we compare the kind of results that can be obtained 
with bounds in the Laplace domain with the results of Sections |3] and [61 In Section [8] we 
state some basic results including bounds on non-homogeneous problems associated to the 
wave equation with different kinds of boundary conditions; these results have been used in 
the previous sections. Finally, Appendix |A] includes the already mentioned treatment of 
some problems related to the wave equation by means of rigorous separation of variables. 

Notation, terminology and background. Given a function of a real variable with 
values in a Banach space X, : M — )■ X, we will say that it is causal when ip{t) = for all 
t < 0. If (y? is a distribution with values in X, we will say that it is causal when the support 
of is contained in [0, oo). The space of fc-times continuously differentiable functions 
/ — >■ X (where / is an interval) will be denoted C'^(/;X). The space of bounded linear 
operators between two Hilbert spaces X and Y is denoted £(X, Y) and endowed with the 
natural operator norm. Standard results on Sobolev spaces will be used thorough. For 
easy reference, see [1] or [H]. Some very basic knowledge on vector- valued distributions on 
the real line will be used: it is essentially limited to concepts like differentiation, support, 
Laplace transform, identification of functions with distributions, etc. All of this can be 
consulted in iTO 1 . 



On time differentiation. There will be two kinds of time derivatives involved in this 
work: for classical strong derivatives with respect to time of functions defined in [0, oo) 
with values on a Banach space X (understanding the derivative as the right derivative at 
t = 0), we will use the notation u; for derivatives of distributions on the real line with 
values in a Banach space X, we will use the notation u'. Partial derivatives with respect 
to t will only make a brief appearance in a formal argument. 

Remark 1.1. If u : [0, oo) X is a continuous function and we define 



then Eu defines a causal X -valued distribution. If u G C^([0, oo);X) and u{0) = 0, then 
[Eu)' = Eu. Also, if u is an X -valued distribution and X C Y with continuous injection, 
then u is a Y -valued distribution and their distributional derivatives are the same, that is, 
when we consider the X-valued distribution u' as a Y -valued distribution, we obtain the 
distributional derivative of the Y -valued distribution u. This fact is actually a particular 
case of the following fact: if u is an X-valued distribution and A G >C(X, Y), then Au is 
a Y -valued distribution and {Au)' = Au' . 





(1) 
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2 Retarded layer potentials 



Let fl- be a bounded open set in M.*^ with Lipschitz boundary T and let := \ 
We assume that the set is connected. No further hypothesis concerning the geometric 
setup will be made in this article. The normal vector field on F, point from Q_ to fl^ will 
be denoted u. 



Claissical integral form of the layer potentials. For densities A, </? : F x R — > M 

that are causal as functions of their real variable (time) , we can define the retarded single 
layer potential by 

and the retarded double layer potential by 

= /^i^fr^ Wy.' - 1== - y|) + 1^ - y|^(y,f - 1^ - y|))dr(y). 

These are valid formulas for x G \ F as long as the densities are smooth enough. The 
two dimensional layer potentials are defined by 

and 

{V.<p)M [ vi^^t-\^-y\) (x-y).Ky) ^^^^ 

STrir |x-y| ^(t - r)^ - |x - yP ' 

1 f r-^^^{y.r) (x-y)..(y) 



27ryrio (t-T)2-|x-y|2 ^(^ _ _ _ y| 



Convolutional notation for potentials and operators will be used throughout. As we will 
shortly see, the convolution symbol makes reference to the time-convolution. 



Layer potentials via transmission problems. In a first step, layer potentials can be 
understood as solutions of transmission problems. Let 7^ (resp. 7+) denote the operator 
that restricts functions on Q_ (resp. to F, i.e., the interior (resp. exterior) trace 
operator. Let similarly denote the interior and exterior normal derivative operators. 
Jumps across F will be denoted 

l^u\ := -f~u - 7+1*, ld^u\ := d^u - d^u, 

while averages will be denoted 

M + ^+u), id^ul := \{d-u + dtu). 
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Given a causal density A, the single layer potential u := S * X can be formally defined as 
the solution to the transmission problem: 



Uu = Au in M'^yr X (0,oo), (2a) 

|7m]=0 onrx(0,oo), (2b) 

|a^M]=A onrx(0,cx)), (2c) 

u{-,0)=0 mR'^\T, (2d) 

Mt(-,0)=0 inM'^\r. (2e) 



Similarly, for a causal density (p, u := V * (p is the solution of the transmission problem: 





Au 


in 


M"^ \ r X (0, oo 


M = 


-if 


on 


r X (o,cx)), 




= 


on 


r X (o,cx)), 


u{ ■ ,0) 


= 


in 




M ■ ,0) 


= 


in 





With this definition, it follows that 

|^(5*A)1 = 0, ld,{V*if)} = 0, |9,(5*A)] = A, t{V*if)j = -if. (3) 

The definition of the layer potentials through transmission problems allows us to define 
the following four retarded boundary integral operators: 



V*A := {{7(5* A)}} =7^(5* A) =7+(5* A), (4) 

/C**A := {R(5*A)}}, (5) 

JC^if := {MV*ip)}}, (6) 

W*^ := -id,{V*f)}} = -d;{V*^) = -d^{V*^). (7) 



These definitions and the jump relations (|3]) prove then that 

d^{S * A) = ^iA + /C* * A 7^(1) *ip) = ±y + IC*ip. 



Layer potentials via their Laplace transforms. Although the definition of the layer 
potentials through the transmission problems they are due to satisfy leads to an easy 
formal introduction of potentials, integral operators and most of the associated Calderon 
calculus with integral operators, properties of these operators are usually obtained by 
studying their Laplace transforms. This is the usual rigorous way of introducing these 
potentials (see |2], [3]). In order to do this, consider the fundamental solution of the 
operator A — for s G C+ := {s G C : Res > 0}: 

|if(')(*s|x-y|), id = 2), 
E,(x,y;s) := { e-^'^-yl 
47r|x — y| 
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The theory of layer potentials for elliptic problems (see |9] or the more general introduction 
in the monograph [H]) can then be invoked in order to define the single and double layer 
potentials, which are weak forms of the integral expressions 

H-'/^{T) s(s) A := Ie,{-, y; ^)A(y) dr(y), 

and 

H'/\r) D{s)^ := ^ VyE,{ ■ , y; s) ■ u{y) ^(y)dr(y), 

respectively. For an arbitrary open set O, we let 

Hi{0) := {u e H\0) : Am G L^{0)], 

endowed with its natural norm. Then S(s) : H-^''^{T) H]^{W^\T) and D(s) : H^''^{T) 
H^{R'^ \ r) are bounded for all s G C_|_. The jump relations 

[7S(.)A] = 0, I9.D(.V1 = 0, |9.S(.)A] = A, tD{s)^} = (8) 

justify the definition of the four associated boundary integral operators using averages of 
the traces 

V(.)A := {{iS{s)X}} = 7^S(s)A, K\s)X := §9.S(s)AS, 

Bounds of the operator norms of the two potentials and four operators above, explicited 
in terms of s, have been obtained in [21 |3] and [12]. Using them, it is then possible 
to use Payley- Wiener's theorem (see an sketch of the theory in [TOj or a full introduc- 
tion in [18]) and show that all six of them (S,D,V, K,K* and W) are Laplace trans- 
forms of operator-valued causal distributions. For instance, it follows that there exists 
an ^(^^"^/^(r), iJ^(M'^ \ r))-valued causal distribution S whose Laplace transform is well 
defined in C+ and is equal to S(s). The theory of vector- valued distributions proves then 
that for any causal if^^/^(r)-valued distribution A, the convolution product iS* A is a well 
defined causal if^(R'^ \ r)-valued distribution. Moreover, if -u := 5 * A, then 

u" = Au. (9) 

(Recall notation for distributional derivatives given at the end of the introductory section.) 
The Laplace operator in is the Laplacian A : HKR"^ \ T) ^ L2(M^ \ T) = L2(R°') and 
([9]) is to be understood as the equality of two L^(M'^)-valued causal distributions. The fact 
that u is causal and that differentiation is understood for distributions defined on the real 
line (as opposed to distributions defined in (0, oo)), encodes the vanishing initial conditions 
( pdj) and (l2ell . The jump properties of S(s) in ([8]) prove then the transmission conditions 
in ([3]). This gives full justification for understanding m = iS * A as a solution of the 
transmission problem ([2]) with time differentiation (and initial conditions) re-understood 
as differentiation of vector-valued distributions. If V and /C* are the causal operator- 
valued distributions whose Laplace transforms are V(s) and K*(s) respectively, then their 
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time convolutions with a given causal density A satisfy the identities (jlj) and ([5]) thus 
identifying the two possible definitions of the time domain integral operators associated 
to the single layer potential. 

The same considerations can be applied for a rigorous definition of the double layer 
potential in the sense of convolutions of vector-valued distributions. Note that both 
layer potentials had been introduced directly (without using the Laplace transform) in 
the three dimensional case in [11], with a theory that cannot be easily extended to the 
two-dimensional case. 

Propagation, uniqueness and integral representation. Finite speed of propaga- 
tion of the waves generated by layer potentials will be a key ingredient in our theoretical 
setting. For our purposes, only the following aspect will be used. Henceforth we take a 
fixed R> such that 

IF C Bo := B{0;R) := {xER'^ : \x\ < R}. (10) 

We also consider the distance between F and dBo: 

5 := min{|x-y| : X G F,y e 95o}. (11) 

For T > 0, we denote Bt '■= B{0]R + T) and we let jt be the trace operator from 
H\Bt\T) to H^/^dBr). 

Proposition 2.1. Let A be an H~^^'^ (T) -valued causal distribution, an H^/^(T) -valued 
causal distribution, and u:=S*\ + T>*ip. 

(a) The temporal support of the H^/'^{dBT) -valued distribution 'Jtu is contained in [T + 
6, oo). 

(b) Letting Or := M'^ \ Bt-5/2, the temporal support of the H^{Ot) -valued distribution 
u\oj. is contained in [T -|- 5/2, oo). 

Proof. This result is a consequence of some simple techniques related to the Laplace 
transform. Firstly, if the Laplace transform F(s) of a distribution / can be bounded as 

||F(s)|| < Cexp(-cRes)|s|^ Vs G C with Res >0, (12) 

where c > and /i G M, then the support of / is contained in [c, oo). Using estimates of 
the fundamental solution Ed as a function of s, it is possible to prove a bound like ( IT2l) 
for S(s) (resp. D(s)) as an operator from H~^/'^{T) (resp. H^/'^{V)) to H^^'^^dBx) and to 
H\Ot). □ 

Proposition 2.2. Let A be an H~^^'^{T) -valued causal distribution and an H^/'^iV)- 
valued causal distribution Lp and assume that both are Laplace transformable. Then u := 
S*\ — 'D*ip is the only causal H^{W^\V) -valued distributional solution of the transmission 
problem 

u" = Au, tu] = if, [d^u] = A 
that admits a Laplace transform. 
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3 Main results 



Before stating the two main results of this paper, we need to make precise statements on 
some constants related to the geometric setting and the Sobolev norms. The reference 
radius R > will be chosen so that ( ITOl) holds. 
Given an open set O, we will denote 



\u\\o ■-- 



o 



N 1/2 / X 1/2 



The first set of constants we need are the Poincare-Friedrichs constants on the balls 
Bt := B{0;R + T) for T > 0: 

\\v\\Br < CrWVvWBr Vt; G H^oiBr). (13) 

A simple scaling argument shows that Ct = Cq{1 + T/R). The second relevant constant 
is a continuity constant for the interior and exterior trace operators. It will be jointly 
expressed for functions that are on each side of F: 

||7^w||i/2,r <C^r||M||i,Bo\r e H\Bo\T). (14) 

Here || • ||i/2,r is a fixed determination of the if ^/^(F)-norm (several equivalent choices 
are available in the literature; see |1] or [I4j|). The third constant is related to a lifting 
of the trace operator. Since 7" : if^(f2_) — )• if^/^(F) is bounded and surjective, there 
exists a bounded linear operator L~ : iJ^/^(F) — H^(n^) such that 'j~L~Lp = ip for all 
Lp E i/^/^(F), i.e., L~ is a bounded right-inverse of the interior trace. We then denote 
Cl := ||i^~||. The use we will make of this operator and its norm will be through L : 

given by 

^= I in Q^, 

noting that 

\\L^p\\i,Rd\r = \\L'p\\i,n_ < CL\\(p\\i/2,r, l~Lip = ip, 7+Lv9 = 0. (15) 

The final constant is related to the definition of the normal derivative. Given u G H\{Bq\ 
F) we can define d'^u G H~^^'^{T) with Green's formula. Then, there is a constant Ci, 
such that 

\\dM-i/2,T < a(||Vw|||„nf,^ + ||A«|||„nn±)'^' V« G Hl{B, \ F). (16) 

The main theorems of this paper are given next. For simplicity of exposition, we 
assume that data are smooth (i.e., C°°) and causal. Some considerations on the smoothness 
of data will be made in Remark 17. 1[ 

Theorem 3.1. Let X be a causal smooth H~^^'^{T) -valued function and let 
B, ''\X,t) := f (||A(r)||_i/2,r+ ||A(r)||_i/2,r)dr. 
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Then for allt>0 

11(5 * A) (t) II 1,1,. < Cr{\\Xm-i/2,r + Vl + c!B-'/\X,t)), (17) 
||(V*A)(t)||i/2,r < C'^{\\Xm-i/2,r + VlTc!B-'/\X,t)), (18) 
||(/C**A)(t)||_i/2,r < v^aCr(||A(t)||_i/2,r + i?2"'^'(A,t)). (19) 
Theorem 3.2. Let Lp be a causal smooth H^^'^iV) -valued function and let 

Bl'\^,t) := £(||¥^(r)||v2,r + ||¥3(r)|K/2,r)dr, 

BT{v,t) := (4||v^(r)||i/2,r + 5||y3(r)|K/2,r + ||¥.(^nr)||iAr)dr. 
Then for allt>0 

||(P*^)(t)||i,i,.\r < CL[Mt)\\,/2,r + V^TcfBl/'{^,t)), (20) 

||(/C*v9)(t)||i/2,r < CrCL{\\ip{t)h/2,r + VT+c!Bl/\cp,t)y (21) 

||(W*y.)(t)||_v2,r < V2C.CL(nvmy2,r + 2\\m\\i/2,r + B'J\^,t)). (22) 



4 The single layer potential 

Since the convolution operator A i— )■ 5 * A preserves causahty, in order to obtain bounds 
at a given value of the time variable t = T, only the value of A in (T, oo) is not relevant. 
Therefore, we can assume without loss of generality that the growth of A allows it to have 
a Laplace transform. We can actually assume that A is compactly supported for the sake 
of the arguments that follow. 



Introduction of a cut-off boundary. Let u := S*X. By Proposition [2]2l m is a causal 
distribution with values in X := H^{M.'^) nif^(M'^\r). Moreover, it is the only (X- valued 
causal distributional Laplace transformable) solution of 

u" = Au and Id^uj = A, (23) 

with the differential equation taking place in the sense of distributions with values in 
L^(M'^ \ r) = L^(M'^), while the transmission condition is to be understood in the sense of 
i?~^/^(r)-valued distributions. Let now T > be fixed and let Bt and S be as in Section 
|2j We look for a causal distribution with values in 

XT:=H^{BT)nHi{BT\T) 

such that 

u't = Aut and IOuUt] = X. (24) 

This differential equation is understandable in the sense of L^(i?T)-valued distributions. 
We will show that for smooth data A, this problem has strong solutions, with the time 
derivatives understood in the classical sense. 
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Proposition 4.1. As H\{Bt \ T)-valued distributions, u = ut in (— oo,T + 5). 

Proof. Consider the H\{Bt \ r)-valued distribution w := u — ut = u\bj. — mt- Then 

w" = Aw, {duw} = 0, and •jtw = Itu. 

Since the support of '-yxu is contained in [T + 5, oo) (by Proposition 12. ip . so is the support 
of w, which proves the result. □ 

Proposition 4.2. For causal A 6 C^(M; if~^/^(r)), the unique solution of fl24p satisfies 
UT e C\[0, oo); L\Bt)) n C\[0, oo); H',{Bt)) n C([0, oo); Xt), (25) 
the strong initial conditions Mt(0) = 'Ur(O) = and the bounds for all t > 

\\uTmi,Br < Cr[\\X{t)U/2,r + ^/l + C^B,'/\X,t)), (26) 
||V«T(t)||B, < Cr{\\X{t)U/2,r + B,'/'{X,t)'), (27) 
WAurmBrW < Cr(||A(t)||_i/2,r + i?2"'/'(A,t)). (28) 

Proof. Consider first the function uq : [0, oo) — )■ Hq{Bt) defined by solving the steady- 
state problems 

-Auoit) + Uo{t) = in Bt\T, ld,Uo{t)} = A(t), 7TMo(t) = 0, 
for t > 0. The variational formulation of this family of boundary value problems is 

" Uo{t) G H^{Bt), 

_ {Vuo{t),Vv)Br + {uo{t),v)Br = (A(t),7^)r Vy G i/o'(^T), 

where (■,-)r is the if~^/^(r) x //^/^(F) duality product. Therefore, a simple argument 
yields 

\\uomi,Br < Cr||A(t)||_i/2,r, || A«o(t) ||B,\r < Cr||A(t)||_i/2,r. (29) 

Note that Uo{t) is the result of applying a bounded linear (time- independent) map H~^^'^{T) — )■ 
Xt to X{t). Therefore, since A is twice continuously differentiable in [0, oo), it follows that 

\\Mt)\\i,Br < Cr||A(t)||-i/2,r, \\Aiio{t)\\Br\r < C^r||A(t)||-i/2,r. (30) 

We next consider the function Vq : [0, oo) — > H'^{Bt) H Hq{Bt) that solves the evolution 
problem 

voit) = Avoit) + uoit) - iioit) t>0, t;o(0) =?)o(0) = 0, (31) 

i.e., the hypotheses of Proposition 18.11 hold with f = uq — uq. Therefore, using (I29l)- (l30l) . 
it follows that 

WAvomB^ < [ \\Vuo{T)-Vuo{r)\\Brdr<CrB^'/\X,t), (32) 
Jo 
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as well as 



VomBr<CTCTB-^'\\t), 



\\VvomBr<CTB-''\\t). 



(33) 



If we now define ut '■= uq + vq, then the regularity requirement ( |25l) is satisfied and the 
three bounds in the statement of the proposition are direct consequences of (l29l) . (l32i) and 
(1331) . Moreover, 

ur(t) = Aurit), Id.UTit)] = X{t), -fTUrit) = Vt > 0. 

Note also that ut{0) = mo(0) = and uriO) = ito(O) = 0, since A(0) = and A(0) = 
(A : R — )■ if~^/^(r) is assumed to be and causal). Therefore, considering the extension 
operator ([T]) it follows that Eut is an X^-valued causal distribution, (Eut)" = Eut = 
EAut = AEut and {d^Eur] = Eld^ur] = ^A|(o,oo) = A. Therefore, Eut satisfies 
and the proof is finished. □ 

Proof of Theorem 13.11 By Proposition 12. the distribution ^Ij^d^ ^^ ^ ^ vanishes in the 
time interval (— oo,T + 6/2). Therefore, by Proposition 14. 1^ UT{t) = in the annular 
domain Bt \ -Bt-<5/2 for all t < T + 6/2. This makes the extension by zero of mt(^) 
toR'^\B^ an element of if^(R'^ \ F) for all t < T + 6/2. (Note that the overlapping 
annular region is needed to ensure that the Laplace operator does not generate a singular 
distribution on dBx.) Then, the argument of Proposition 14.11 can be used to show that 
the distribution u can be identified with this extension in the time interval (— oo, T + 5/2). 
Therefore, identifying u{T) = ut{T), the inequalities of Proposition 14.21 yield 



We can now substitute all occurrences of T by t, since T was arbitrary. The result is now 
almost straightforward. First of all, is just f lT7|) . Also, by the trace inequality f lT^ 
and the fact that V * A = 7^(5 * A), (HH]) is a direct consequence of f lT7|) . Finally, the 
bound for the normal derivative (fT6|) . the fact that /C* * A = ^du{S * A)^, and inequalities 
(I35])-(|36D prove 

5 The double layer potential 

We start by introducing a cut-off boundary OBt as in Section H] (for arbitrary T > 0). 
We are going to compare u := V * ip with the causal distribution ut with values in 




(34) 



(35) 



(36) 



YT:={veHi{BT\T) : JtU = 0}, 



such that 



u'^ = Aut, Ijut} = ~f and 



IduUr] = 0. 



(37) 
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The same argument as the one of Proposition 14.11 shows that, as H\{Bt \ r)-valued 
distributions u = ut in (— oo,T + 5), where 5 is defined in ( ITTj) . Smoothness of the 
solution of (1371) and bounds for it in different norms will be proved in two steps. Note 
that, from the point of view of regularity Proposition 15.21 improves the initial estimate of 
Proposition 15.11 but that more regularity of if is used in the process. 

Proposition 5.1. For causal ip G C^(M; if^/^(r)), the unique solution of (l37l) satisfies 

UT G Ci([0, oo); L\Bt)) H C([0, oo); H\Bt \ P)), (38) 
the strong initial conditions mt(0) = uri'^) = and the hounds for all t > 

\\uTmi,Br\r < CL{Mmi/2,r + + Bl^'{v,t)), (39) 
||VwT(t)b.\r < CL{\Mt)\\,/,,r + Bl^'icp,t)y (40) 

Proof. Let first uq '■ [0, oo) — )• H^{Bt \ P) be given by solving the steady-state problems 

-Auo{t) + uo{t) = in 5t \ P, buoit)] = -if{t), (41a) 
7TMo(t) = 0, ld,uo{t)j = 0, (41b) 

for each t > 0. The variational formulation of ( 14T1) is 
- Uo{t)eH\BT\T), 

tuoit)] = -ifit), 7T«o(t) = 0, (42) 
_ (VMo(t), ^v)Br\r + iuoit),v)Br =0 \fv E HI{Bt). 

Using the lifting operator (ITS!) . we can choose the test v = Uo(t) + Lip{t) G Hq{Bt) in 
(H2|) . and prove the estimate 

\\uo{t)h,Br\r < \\L^mi,Br\r < CL||^(t)||iAr. (43) 

Since uo{t) is the result of applying a hnear bounded (time- independent) map if^/^(P) — )■ 
Yt to (p{t), it follows that 

\\uomi,Br\r < CL\\m\\i/2,r. (44) 

We then consider vq : [0, oo) — i- Hq{Bt) to be a solution of 

Vo{t) = Avoit) + uo{t) - iioit) t>0, vo{0) = Vo{0) = 0, (45) 

with the equation taking place in H^^{Bx) (that is, Vq is a weak solution in the terminology 
of Section [8]). By Proposition 18.21 (the right-hand side f := uq — uq : [0, oo) — ?■ L'^{Bt) is 
continuous) we can bound 

\\Vvo{t)\\B,< [ \\uo{T)-iio{r)\\Brdr<CLBl/\v,t), (46) 
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where we have apphed ( 143|) -( 14^ . 

Let us then define ut '■= uq + vq. Since (p G C^([0, oo); H^^'^{T)), it follows that uq G 
C^{[0, oo) ■,H\Bt\T)) and vq G C\[0, oo); L^Bt)) D C{[0, oo); H^{Bt)) by Proposition 
18.21 Therefore, ut satisfies fl55]l . Since ip{0) = 0(0) = 0, it follows that mt(0) = 'Ur(O) = 0. 
Considering fHTl) and fH51) (recall that fo takes values in Hq^Bt)), it follows that 

UT{t) = ^Mt), buTm = -^{t), ld,UTm = 0, 7t«t = 0, Vt>0. 

Noting that ||^^o(^)||bt < C'tII Vt;o(t) Hbj,, and using (gS]), dS]), and it follows that 
satisfies the bounds fl39l) and fj40l) . 

The delicate point of this proof lies in showing that ut can be identified with the 
Ir-valued distributional solution of ( l37l) . since vq is not a continuous 1^- valued function. 
However, wo(t) := fo(T)dr is a continuous function with values in H'^{Bt) n Hq{Bt) 
(see Proposition 18. 2p and therefore in Yt. We can then define ut '■= Euq + (Ewo)', which 
is a causal Ir-valued distribution for which we can easily prove that 

tur] = Etuo] + (EtwolY = -^<^|(o,oo) = 

and similarly [d^UTj = 0. Since Wq G {[0, oo); L'^{Bt)) n C([0, oo); Ft), and Wo(0) = 
Wo{0) = 0, it follows that (Ewq)" = Ewq = EAwq = AEwq and therefore u'^ = (Euq)" + 
(EAwq)' = AEuq + A{Ewo)' = Aut- Thus, ut satisfies (jnTj). Finally, since Wq G 
{[0, oo); Hq{Bt)) and Wo{0) = 0, it is clear that, as an ifQ(_B7^)-valued distribution 
{EwqY = Ewq = Evo and thus, as an H^[Bt \ r)-valued distribution ut = Eut, and the 
bounds fl39|) and fHOl) are satisfied by the solution of fl37|) . □ 

Proposition 5.2. For causal Lp G C^{E^; H^^'^iV)), the unique solution of ( 1371) satisfies 

UT G C2([0, oo); L\Bt)) H C1([0, oo); //^(Et \ F)) n C([0, oo); Yt)) (47) 
and the hounds for all t > 

||AnT(t)b,\r<^^L(4||v5(t)l|i/2,r + 2||0(t)||i/2,r + i?r(^,t)). (48) 
Proof. Consider now the solution of the problems 

-Aui{t) + ui{t) = L{Cp{t) - ip{t)) in Bt \ F, hui{t)\ = -ip{t), (49a) 

7TWi(t) = 0, ld,uiit)j = 0, (49b) 

for each t > 0, where L is the lifting operator of f[T51) . Using the variational formulation 
of dH]) and the fact that ui{t) + Lip{t) G H^{Bt), it follows that 

||^ii(t)+L^(t)||i,B,\r < \\L^{t)h,B^\r+\\L{m~Vm\Br < CL(2||<^(t)||i/2,r+||</3(t)||i/2,r) 
and therefore 

\\Mt)hBr\r < CL(3||v^(t)||i/2,r + WmWmr). (50) 
Using (iH]) and ([50]) it also follows that 

||Ani(t)||B, < \Mt)\\Br + mm - m)\\Br < C^l(41| ^(t) || i/2,r + 2||0(t)|K/2,r). (51) 
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Differentiating ( ISOj) twice with respect to t, it follows that 

\\Mt)\\i,Br\r < CLmm\\i/2,r + y\t)\U/2,r). (52) 
Consider next the evolution equation that looks for vi : [0, oo) — )■ Yt such that 

vi = Avi{t) + f{t) Vt>0, t;i(0) = t;i(0) = 0, (53) 

where 

fit) := ni(t) - ni(t) + L(v9(t) - v3(t)) = Ani(t) - ni(t). 

Note that |7/(t)] = for all t, and that / : [0, oo) — > Hq{Bt) is continuous. Moreover, 
by (l50ll and (1521) . we can bound 

W^mU < ll/(t)l|i,B, <C,(4||v,(t)||v2,r + 5||y3(t)|K/2,r + l|^(^nt)lli/2,r). (54) 
By Proposition I8.H problem fl53|) has a unique (strong) solution and we can bound 

\\^v,{t)\\B, < f ||V/(r)b,dr < CLB,{^,t). (55) 

If we finally define ut '■= wi + fi, the smoothness of ui : [0, oo) — )■ Ft (directly inherited 
from that of if) and the regularity of vi that is derived from Proposition 18.11 prove that 
P7|) holds. The bound psj) is a direct consequence of flSTl) and fl53|) . The fact that the 
extension Eut is the Ir-valued causal distributional solution of fl37|) can be proved with 
the same kind of arguments that were used at the end of the proof of Propositions 14.21 □ 

Proof of Theorem 13. 2L With exactly the same arguments that allowed to prove ( l34l) . 
(1351) and ( l36l) as a consequence of Proposition 14.21 we can prove that for allt > 

\\{V * ^){t)\\,^^.\r < CL{\W)\\il2.v + V^+ClBl'\^,t)), (56) 
||V(I?*^)(t)|U.\r < CL{\\v{t)\\y2.y + Bl'\v,t)), (57) 
||A(P*^)(t)|U.\r < CL[AMt)\\i/2,v + 2m)\\i/2,v + BT{^.t)). (58) 

as a consequence of Propositions 15.11 and 15.21 The bounds of Theorem 13.21 are now 
straightforward. Inequality fl20l) is just fl56|) . while the fact that }C* ip = ^'y(T> * ip)^ and 
the trace inequality f[T^ prove fl2U]) . Finally the bound for the normal derivative flTBl) . 
the definition of W * = — (9;^(P * (/?) and inequalities flS71)-flS51) prove fl2^ . 

6 Exterior Steklov-Poincare operators 

In this section we include bounds on the exterior Dirichlet-to-Neumann and Neumann-to- 
Dirichlet operators that can be obtained with the same techniques than in the previous sec- 
tions. We give some details for the easier case (the Neumann-to-Dirichlet operator,whose 
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treatment runs in parallel to that of the single layer retarded potential) in order to em- 
phasize the need of dealing with some slightly different evolution problems as part of the 
analysis process. 

Let us consider the bounded open set B:^ := Bt H Q+ and the space 

Vt := {u G H\B+) : = 0}. (59) 
We can then consider the associated Poincare-Friedrichs inequality 

\\u\\s+ < Et\\Vu\\^+ Vm e Vt. (60) 

Recalling that Bt is a ball with radius R + T, it is possible to take Et < 2{R + T) (see [71 
Chapter II, Section 1]). Since the exterior trace operator 7"*" : Vq — )■ H^^'^{T) is surjective, 
it has a bounded right-inverse. By extending this right-inverse by zero to ^2+ \ Bq we can 
construct L+ : H^/^{T) H^{n+) satisfying 

||L+(^||i,n+ = ||I^+V^|li,B+ < C+||<^||i/2,r, l^L+y^ = ^ V<^ G H'/\r). (61) 

Note that, in particular, jtL'^'^ = for all T > and all (p. 

Theorem 6.1. For causal A e C^{R; H-^/^(T)), the unique causal H^{Q^)-valued Laplace 
transformable distribution such that 

u" = Au, dtu = A, (62) 

satisfies the bounds 

Mt)\\i,n+ < Cr(||A(t)||_i/2T + ^Jl + E^B^'/\X,t)^ Vt > 0. (63) 

Finally the associated Neumann-to-Dirichlet operator NtD(A) := j'^u (where u is the 
solution of ( l62i) ) satisfies the bounds 

||NtD(A)(t)||i/2,r < Cl (11 A(t) || _i/2,r + sJl + El B'^'^X, t)) Vt > 0. (64) 

Proof. The proof is very similar to the one of Theorem 13.11 By solving steady state 
problems, we first construct uq : [0, 00) — )■ H^{B^) satisfying 

- Auo{t) + uo{t) = in 5+, d^uo{t) = A(t), jTUo{t) = 0. (65) 

A simple argument allows us to bound 

ho(t)||i,B+ < C^r||A(t)||_i/2,r and ||«o(t)||i,B+ < Cr||A(t)||_v2,r. (66) 

This function feeds the evolution equation looking for vq : [0, 00) — )• Dt (the set Dt is 
defined in (|79ll ) that satisfies 

vo{t) = Avoit) + uoit) - iio{t) t > 0, vo{t) = vo{t) = 0. (67) 
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We now apply the general result on the wave equation with mixed boundary conditions 
(Proposition I8.3p that guarantees the existence of a strong solution of (!67|) satisfying the 
bounds 

\\vo{t)h^ < ETCrB-'/\X,t), WVvoim^^ < CrB~''\\t). (68) 

Adding the solutions of (1651) and (1671) we obtain a function ut '■= Uq + : [0, oo) — i- 
H\{B^) n Vt satisfying u{t) = Au{t), dyUrit) = A(t) and vanishing initial conditions 
at t = 0. The extension Eut is then an {H\{B^) fl VT)-valued causal distributional 
solution of ([S3) (with the Laplace operator acting only in the bounded domain B^). 
The arguments of Proposition 14.11 and at the beginning of the proof of Theorem 13.11 can 
be applied verbatim in order to identify the function that extends urit) by zero to the 
exterior of Bt with the iJ^(fi+)- valued distributional solution of fl62l) for t < T + 5/2. 
Finally, the bound fl63|) for t = T is a straightforward consequence of fl66l) . fl68l) . and 
the identification of u{T) = ut(T), while ( !64l) follows from ( l63l) and the trace inequality 

m- □ 

Theorem 6.2. For causal G C^(M; H^/^(T)), the unique causal H\{yt^) -valued Laplace 
transformable distribution such that 

u" = Au, = ip, (69) 

satisfies the bounds 

IHt) < Ci^Mt) II i/2,r + ^l + El, Bl'\^, t)) Vt>0. (70) 

Finally the associated Dirichlet-to- Neumann operator DtN(y9) := d:^u (where u is the 
solution of ( 1^ ) satisfies the bounds 

||DtN(^)(t)||_i/2,r< v^aC^^(4||v.(t)||i/2,r + 2||y3(t)||i/2,r + i?r(^,t)). (71) 



Proof. This result can be proved like Theorem 13. 2l by resorting to a double decomposition 
of a localized version of problem ( l69i) (obtained by adding a boundary condition 7tm = 0) 
as a sum of an adequate steady-state lifting of the Dirichlet data plus the solution of 
an evolution problem. The proof is almost identical to that of Theorem 13.21 the main 
difference is in the evolution problem, that now contains Dirichlet boundary conditions 
on r as well as on OBt (see Remark 18.11) . □ 



7 Comparison with Laplace domain bounds 

The original analysis for the layer operators and associated integral equations, given in 
[2] and [3], was entirely developed in the resolvent set (that is, by taking the Laplace 
transform). Those results can be used to derive uniform bounds similar to those of 
Theorems 13.11 and 13. 2[ We next show how to obtain these estimates and show that our 
technique produces stronger estimates, in terms of requiring less regularity of the densities 
and having constants that increase less fast with respect to t. 
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Remark 7.1. Before moving on, let us emphasize that Theorems \3.1\ and \3.2\ do not 
require the highest order derivative involved to be bounded and local integrability is enough 
to keep the quantities 5^^^^ and bI^"^ bounded. This is as much as saying that Theorem 
I3l\ IS still valid if \ : R ^ H-^/^{T) is causal and (therefore A(0) = A(0) = 0) and 
X" is locally integrable. Similarly, estimates fl20p and (1211) of Theorem IJ.jjl hold for 
causal densities with locally integrable second derivative, while the estimate (1221) holds for 
densities with locally integrable fourth derivative. 

Estimates in the Laplace domain can be obtained using the following all purpose 
theorem, which is just a refinement of Lemma 2.2 in [13]. The refinement stems from 
taking more restrictive hypotheses; these ones are chosen in order to fit closer to what 
can be proved for all operators associated to the wave equation. 

Theorem 7.1. Let f be an C{X,Y)-valued causal distribution whose Laplace transform 
F(s) exists for all s G C+ := {s G C : Res > 0} and satisfies 

\\F{s)\\^x,Y) < CA^^s)\sr VsgC+, 

where /i > and C-p : (0, oo) — )■ (0, oo) is a non-increasing function. Let 

k:=[fx + 2\, e:=k-{fx + l) e{0,l]. 

Then for all causal C^^^ function g : ^ X with locally integrable k-th distributional 
derivative, the Y -valued distribution f * g is a causal continuous function such that 

Jo 

In Table [7] we compare regularity and growth of the bounds between what Theorems 
13.11 and 13.21 prove and what can be obtained by a systematic analysis in the Laplace 
domain. The bounds in the Laplace domain are explicit or implicitly given in [2] and 
[3J. They are also collected in Appendix 2]. We also include the comparison of 
what Theorems 16.11 and 16.21 assert about Steklov-Poincare operators with similar results 
obtained through the Laplace domain analysis. 



8 Basic results on some evolution equations 
8.1 Homogeneous Dirichlet conditions 

In this section we gather some results concerning solutions of the non-homogeneous wave 
equation with homogeneous initial conditions and homogeneous Dirichlet boundary condi- 
tions on the ball B^ introduced in Section [2j We recall that Ct is the Poincare-Friedrichs 
constant in Bt (see ( fT3l) ). The problem under consideration is: 

u{t) = Au{t) + f{t) t > 0, (72a) 
^Tu{t) = t>0, (72b) 

u{0) = m(0) = 0. (72c) 
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X 



CF(a)|s|^ E{t)D,{rj,t) 



n 0{t) 



S H-^'\T) ^ H\W^) 

K* iy-V2(r) ^ i7-i/2(r) 

D H^/\V) ^ H\W^\V) 

w iyi/2(r) ^ iJ-V2(r) 

NtD H-^/^iV) H^/^iV) 

DtN H^I^V) ^ i/-V2(r) 



(JC 



3/2 



<T<£372 
<To;3/2 

(To:3/2 



max{l,t2}/}3(A,t) 
max{l,t2}/}3(A,t) 
t3/2 max{l,t3/2}/)3(A,t) 

t3/2 max{l,t^/^}D3{ip,t) 
t3/2 max{l,t3/2}/}3(<^,t) 

max{l, t}D4{(p, t) 
max{l,t2}L)3(A,t) 
max{l, t}D4(V9, t) 



0(1) 

0{t) 
0{l) 
0{t) 
0{1) 



Table 1: Bounds obtained using Theorem l7.1l and known estimates in the Laplace domain. 
The first line acts as a prototype and has to be read as ||F(s) ||£(x,y) < C x Cf{(7)\s\^ and 
IK/ * v){t)\\Y <C X E{t) X /J ||r/(^)(r)||xdr. Here a := Res and q_ := min{l,a}.The last 
two columns contain information given by Theorems 13.11 13. 2[ 16.11 and 16. 2[ indicating the 
highest order n of differentiation of t] involved in the bounds for f * r] and the growth of 
the bound clS cl function of t for large t. 



We will deal with two different types of solutions of this problem. A strong solution is a 
function such that 

u G C2([0, oo); L\Bt)) H ^^([0, oo); HI{Bt)) H C([0, oo); H\Bt)). (73) 

with the wave equation satisfied in L'^{Bt) for all t. We will refer to a weak solution as a 
function such that 

u e C2([0, oo); H-\Bt)) n Ci([0, oo); L\Bt)) n C([0, oo); HI{Bt)) (74) 

and such that the wave equation is satisfied in H~^{Bt) (the dual space of Hq{Bt)) for 
all t. Note that the concept of weak solution relaxes both time and space regularity 
requirements and does not exactly coincide with the concept of mild solution given in [16] 
for example. 

Proposition 8.1. Let f : [0, oo) — t- Hq{Bt) be continuous. Then, problem f l72|) has a 
unique strong solution satisfying the bounds for all t > 

\Ht)\\Br<CT f \\fiT)\\B,dT, (75) 

Jo 

||V«(t)||B, < / ||/(r)||B,dr, (76) 
\\Au{t)\\B,< f\\Vfir)\\B,dT. (77) 

^0 
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Proposition 8.2. Let f : [0, oo) — L'^[Bt) be continuous. Then problem ( I72|) has a 
unique weak solution, and the bound (!76|) is still valid. Finally, the function w{t) := 
Jq M(r)dr is continuous from [0, oo) to H'^{Bt). 

Remark 8.1. Propositions \ 8. 1\ and \8.2\ still hold for the Dirichlet problem in the domain 
:= Bt n with the following modifications: the space H'^{Bt) has to be substituted 
by H\{B^), and the constant Ct in fl75|) has to be substituted by the constant Et of the 
Poincare-Friedrichs inequality ( l6Ql) . 

8.2 Mixed conditions 

Let us now consider the set B^ := B^ H Vt^ and the evolution problem 



u{t) = Au{t) + f{t) t > 0, (78a) 

^Tu{t) = t > 0, (78b) 

dtu{t) = t > 0, (78c) 

u(0) = u(0) = 0. (78d) 

We thus consider the spaces Vt given in (|59|1 and 

Dt := {ueVr : Aue L'^{B+), d+u = 0} (79) 



= {uEVtH Hl{B+) : (Vm, Vv)u+ + (m, v)n+ = Vt; G Vt}. 

Proposition 8.3. For f G C([0, oo); V^), the initial value problem fl78|) has a unique 
solution 

u G C2([0, oo); L2(r])) n C^([0, oo); Vt) H C([0, oo); Dt), 

satisfying 

\\u(t)\\B^<CT j\\f{r)\\^^^dT, (80) 
l|V^t)L.<£||/(r)||^+dr, (81) 
||A^t)||^+<^*||V/(r)||^+dr. (82) 

If f & C([0, oo); L^(fi)) there exists a unique weak solution of ( ITHl) (that is, with the 
equation satisfied inV^) 

u G C2([0, oo); 14) n Ci([0, oo); /.^(fi)) n C([0, oo); I^t), 

satisfying flHTl) anc? swc/i t/iat iy(t) := M(r)dr zs m C([0, oo); Dy). 
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A Wave equations by separation of variables 



In this section we are going to give a direct proof of a generalization Propositions 18.11 
and 18. 2[ This proof will be based on direct arguments with generalized Fourier series 
and will allows us to obtain the needed uniform estimates of non-homogeneous evolution 
equation of the second order in terms of norms of the data. The Hilbert structure of 
the functional spaces is going to be used in depth, allowing us to obtain strong results 
that cannot be easily derived with a direct application of the best known results on the 
theory of semigroups of operators. This is not to say that these results do no exist, but 
we think it can be of interest (especially within the boundary integral community) to see 
a direct proof of these theorems based on functional analysis tools that are common for 
researchers integral equations. 

A.l Three lemmas about series 

In all the following results X is a separable Hilbert space and / := [a, b] is a compact 
interval. 

Lemma A.l. Assume that Cn '■ I ^ X are continuous, 

(c„(t),c„(t))x = Vn^m, Vt G /, (83) 

and 

\\cn{t)\\x < Mn \/teI, Vn With ^M„<CX). 

n=l 

Then the series 

oo 

c(t):=5^c„(t) (84) 

n=l 

converges uniformly in t to a continuous function. 
Proof Let Sn := En=i Cn e C(/; X). For all M > N, 

M M 
\\sM{t)-Sr,{t)\\\= l|Cn(t)||^< Yl 

n=N+l n=N+l 

which proves that SN{t) converges uniformly. Continuity of the limit is a direct conse- 
quence of the uniform convergence of the series. □ 

Lemma A. 2. Assume that Cn '■ I —> X are continuously differentiahle, 

(c„(t),c^(r))x = Mn^m, yt,T e I, (85) 

and 

oo 

\\Cnit)\\x + \\Cn{t)\\x < Mn ^t & I , Vn With ^ M„ < CX). 

n=l 

Then the uniformly convergent series (l84l) defines a C^(/;X) function and it can be dif- 
ferentiated term by term. 
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Proof. The hypothesis (!85|) imphes ( 183|) as well as 

(c„(t), Cmit))x = Vn ^ m, Vt G /. 
If SN:=En=iCneC\l;X),then 



\\sM{t)-SNmx + \\sM{t)-SNmx= Yl l|CnWllx< I] M 



M M 

A// 



n=N+l n=N+l 

and therefore sat is Cauchy in C^{I; X) and thus convergent. The fact that the derivatives 
of the series converges to the series of the derivatives is part of what convergence in 
J; X) means. □ 

Lemma A. 3. Let f : I ^ X be a continuous function and let {0n} be a Hilbert basis of 
X. Then 



f{t)=J2ifit),<Pn)xC 



n=l 



uniformly intEl. 



Proof. Note first that for fixed t, f{t) G X can be expanded in the Hilbert basis, so 
convergence of the series is easy to prove. Next, consider the square of the norms of the 
A^-th partial sums 



N 2 ^ 



aN 



1 1 X 



which are continuous functions of t. The pointwise limit is which is also a contin- 

uous function of t. Since the sequence is increasing, by Dini's Theorem, convergence 
o-N ||/( ■ is uniform. Finally 



^ 2 



fit) -Y^Uit). <Pn)x<Pn ^= Y |(/(t),0n)xr = ||/(t)||^-a;v(t), 
n=l n=N+l 

which proves the uniform convergence of the series. □ 

A. 2 The Dirichlet spectral series of the Laplace operator 

Let f2 be a Lipschitz domain and consider the sequence of Dirichlet eigenvalues and 
eigenf unctions of the Laplace operator: 

0„ G Hq{Q) - A0„ = A„0„. 

The sequence is taken with non- decreasing values of A„ and assuming (f)m)n = ^nm, 
for all m, n, i.e., L^(n)-orthonormality of eigenfunctions. Thus, {</)„} is a Hilbert basis of 
L^(fi) and consequently, for all u G L^(fi) 



Ml 



$^IK0n)nr (86) 



n=l 
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and 

oo 

U = ^{u,(f)n)n(pn, (87) 

n=l 

with convergence in L^(fi). Using the orthogonahty (V^^, V(f)rn)n = ^nm^n, we can prove 
that 



71=1 

and 

oo 



ra=l 



This expression gives a direct estimate of the corresponding Poincare-Friedrichs inequahty 

as 

ll^llf^ < , ^ ^ IIVmII^ =: Co||V«||q Vm G i/o^(f^). 
Vmm An 

Moreover, if m G Hq{Q), the series representation ( l86i) converges in iJQ(fi). 
The associated Green operator is the operator G : L'^{Q) — > D{A) given by 

u := Gf solution of u G Hq{^1), —An = / in fi. 

Here -D(A) := {u G -f^o(^^) : Am G L^(fi)}. Note that for the case of a smooth domain 
D{A) = Hq{Q) n ahhough this fact will not be used in the sequel. The space 

D{A) is endowed with the norm ||A ■ The series representation of G is given by the 
expression 



Gf = J2K\f,' 



n=l 



(with convergence in ((])). Picard's Criterion can then be used to show that G is 
surjective and 

oo 

D{A) = {ueL\Q) : 5^A^|(«,0„)np}. 

n=l 

Two more series representations are then directly available, one for the Laplacian 

oo 

-Am = ^ A„(u, Wu g D{A), 

n=l 

(with convergence in L'^{Q)) and another one for its norm 

oo 

l|A«il^ = Y,Xl\{u,<Pnh\' Vm G D{A). (89) 



n=l 
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A. 3 Strong solutions of the wave equation 

We start the section with a reminder of one of the possible versions of Duhamel's principle 
that will be useful in the sequel. Its proof is straightforward. 

Lemma A. 4. Let g : [0, oo) W be a continuous function, u > and define 

a{t) := / u^^ sm{u{t — T))g{r)dT. 
Jo 

Then a G C^{[0, oo), a{0) = d(0) = 0, 

a{t) = / cos(a;(t — T))g{r)dT 
Jo 

and a(t) + u'^a{t) = g(t) for all t >0. 

For notational convenience, we will write ■= V^n- 

Proposition A. 5. Let f : [0, oo) — i- Hq{^) be a continuous function and consider the 
sequence 

Unit) := { I E:^^ sin i£„it - t)]( fir), (i)r,)odT ] (t>„, n>l. 
Then, the function 

satisfies 



Tc'sin (a(t-r))(/(r),0„)cdr 
Jo 



u{t) := ^u„(t) 



n=l 



u e C\[0, oo); L\n)) n C\[0, oo); H^{n)) n C([0, oo); D{A)). 
Moreover, u is the unique strong solution of the following evolution equation: 

u{t) = Au{t) + fit) Vt > 0, m(0) = m(0) = 0. 



(90) 
(91) 
(92) 



Proof. As a direct consequence of Lemma [A. 41 it follows that u„ G C^([0, oo); X), where 
X is any of ^^(fi), H^i^) or D(A). Also, for alH > 



Unit) 



COS 



{Ut-r))ifir),(j)n)ndr)(j)n 



Unit) = if it), 0n)c0n " KUnit) = if it), 0n)f70n + AM„(t), 



and 



(u„(t), Mm(r))f7 = (VM„,(t), Vumir))n = 7^ m, Vt, r > 0. (93) 

By dHl]), it follows that for t G [0,T], 

\\^Mt)\\l = / C'sin(a(t-r))(/(r),</.„)^dr 
Jo 

< Xnt f \ifiT),<Pn)n\'dr<T[ A„|(/(r), 0„)f,|2dr =: M«. 
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By the Monotone Convergence Theorem and (l88ll . we easily show that 

oo „j' oo „J' 

Y.^n^=T {Y.>^n\{f{T)AM'yr = T ||V/(r)||^dr. 

n=l -^0 ^ Jo 

Thanks to these bounds and fl93l) (recall that Aun{t) = — A„u„(t)), Lemma [A. II can be 
now applied in the space X = D{A) and interval / = [0,T] for arbitrary T > and we 
thus prove that u G C([0, cx)); L'(A)) C C{[0, oo); H^{n)) C C([0, oo); L2(fi)). Note that 
the series ( l90l) converges for all t and therefore, using the fact that Un{0) = 0, it follows 
that u{0) = 0. Note also that in particular 



Au{t) = J2^Mt) 



(94) 



n=l 



uniformly in t G [0, T] for arbitrary T. 
In a second step, we use flSSl) to bound 



2 +||VMn(t)||Q = A„ 

+A„ 

< T 



fc'sin (a(t-T))(/(r),0„)ndr 
cos {Uit - r)) (/(r), 0n)ndr 

l + A„)|(/(r),0„)d'dr=:Mf . 

By the Monotone Convergence Theorem and the series representations of the norms (186!) 
and (188|) . we obtain 



2) 



T 



n=l 



l|V/(r)||^ + ||/(r)||Mdr. 



Using ( l93ll . we can apply Lemma [A. 2 1 in the space X = -ffo (fi) and the intervals / = [0, T] 
to prove that u G C^([0, oo); El{^)) C C^([0, oo); ^^(fi)). From this, it follows that m(0) = 
0. 

In a third step, we notice that by (IMl) and Lemma lA. 31 



J] (At.„(t) + (/(t) 



71=1 



Aw(t) + /(t), 



with convergence in L^(f2) uniformly in t G [0,T] for any T. Since 'iin(^) = Aunif) + 
(/(t), 0„)o0ri, it follows that the series of the second derivatives L^(f2)-converges, uni- 
formly in t, to a continuous function. Since the series of the first derivatives is t-uniformly 
L^(f2)-convergent (it is actually ifg (^)"Convergent, as we have seen before), it follows that 
u{t) = Au{t) + f{t) for all t > 0, and that u G C([0, oo); L'^{n)). 
Finally, if u satisfies (lOTll and the homogeneous wave equation 

u{t) = Au{t) Vt > 0, u{Q) = u{0) = 0, (95) 

then, a simple well-known energy argument shows that u = 0, which proves uniqueness 
of strong solution to (J92|l . □ 
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Proposition A. 6. Let u he the function of Proposition [A7^ Then, for all t > 0, 

||A«(t)lb< /V/Wlbdr and ||Vn(t)|b < f \\f{r)\\^dT. (96) 



Proof. For arbitrary t > consider the functions ■ ; t) : [0,t] — > -D(A) given by 

^n(T; t) := sin(^„(t - r))(/(r), 

These functions are mutually orthogonal in -D(A) and HI{Vl). Note that V'n := A^^^^^n 
is a complete orthonormal set in HI{VL) and that 

It is then easy to prove the bounds 

||A^?n(r;t)||^ < |(V/(r),VV^„)f^r VrG[0,t], Vn, (97) 
||V^?„(r;t)||^ < |(/(r),0Of,|2 VrG[0,t], Vn. (98) 

Note that by Lemma [A. 31 the series 

oo oo 

5^|(V/(r),V^.)d' = l|V/(r)||^, and |(/(r), 0„)d' = (99) 

n=l n=l 

converge uniformly in r G [0,t]. Using (!97|) and (!99|) . it is clear that 

oo 

[0, t]3T^ g{r- 1) := 9n{r; t) (100) 

n=l 

is well defined as a /)( A) -convergent series. Since convergence of the series (p9l) it also 
follows that the series (llOOl) is r- uniformly convergent in D{A) and therefore in Hq{Q) 
as well. Uniform convergence then allows to interchange summation and integral signs in 
the following equalities 

u{t) = ^Un{t) = ^ / 5'n(^;i)dr= / ^5(„(r;t)dr= / 5((r;t)dr. 

n=l n,=W0 -^0 n=l 

Applying now ( 197|) . ( I99i) . and Bochner's Theorem in the space -D(A), it follows that 

||A^t)|b< r||A^?(r;t)|bdr< T || V/(r) l^dr. 
Jo Jo 

Similarly, (l98l) . (l99l) . and Bochner's Theorem in HliVt)^ prove that 

\\Vu{t)U< f \\Vg{T;t)\\ndT < f \\f{r)\\ndT, 
Jo Jo 

which finishes the proof. □ 
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A. 4 Weak solutions of the wave equation 

In this section we deal with solutions of the evolution problem ( p2|) when / : [0, oo) — )■ 
L^(f2) is continuous. In this case, we will understand the wave equation as taking place 
in H'^{Q) for all t > 0. We first make some precisions about dual spaces and operators. 

As customary in the literature, we let H^^{Q) be the representation of the dual space 
of Hq{Q) that is obtained when L'^{Q) is identified with its own dual space. If we denote 
by ( ■ , ■ )q the corresponding representation of the H~^{Q) x Hq{Q) duality product as 
an extension of the L^(f2) inner product, then 

II II {v,u)n ,_i|^ ^ ^ |2^^/2 

\\v\\^i := sup — ' ^ ^ ^ I ' 



{Y.^n'\{vAnh\') . (101) 



n=l 



The Laplace operator admits a unique extension A : Hq{Q) — t- H ^(i7) given by the 
duality product 

-(Am, v)n = (Vm, Vv)n Vm, v G H^{n) 
and admitting the series representation 

oo 
n=l 

with convergence in H^^{Q). Here A is just the distributional Laplace operator. 

Proposition A. 7. Let f : [0, oo) — )• L^(f2) be continuous. Then the initial value problem 
has a unique solution with regularity 



u e C2([0, oo); H-\n)) n C\[0, oo); L^{n)) n C([0, oo); H^{n)). (102) 
This solution satisfies the bound 



Vu{t)\\n< [ ||/(r)|bdr Vt > 0. 
Jo 



Finally the function w{t) := J^u{T)dT is continuous from [0, oo) to -D(A). 
Proof. Consider first the operator 

oo 

G'/'f:=J2X-'/\fAnh<Pn. 

n=l 

Because of the series representation of the norms (see fl86|) . fl88|) . fl89|) and fllOll) ) it is 
simple to see that G^/^ defines an isometric isomorphism from H~^(Q) to L^(r2), from 
L^{n) to H^{n) and from H^{n) to D{A). It is also clear that AG-^/^ ^ Q-i/^^ ^s a 
bounded operator from -D(A) to H^^{Q). 

As a simple consequence of the above, G^^'^f E C([0, oo); Hq(^1)) and the problem 

i)(t) = Av{t) + G^/V(i) Vt > 0, v{0) = v{0) = 
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has a unique strong solution by Proposition IA.5I We next define u := G^^^'^v. By the 
relations between the norms given by G^^^^ and by the regularity of v given by Proposition 
IA.5t it follows that u satisfies (I102p . It is also clear that u{0) = u{0) = 0. Additionally, 

u{t) = G-^/%{t) = G-'^\Av{t)+G^^^f{t)) = AG-^/\{t)+f{t) = Au{t)+f{t) Vt > 0, 

which makes u a weak solution of ( l92l) . Also, by Proposition IA.61 

\\Vu{t)\\n = \\Avit)\\n< T || VG^/V(r)|bdr = f \\f{T)\\ndT. 

Jo Jo 

To prove uniqueness of weak solution, we note that if u satisfies fll02p and the initial 
value problem fl95|) (with the equation satisfies in H~^{il)), then G^/'^u is a strong solution 
of fl95l) and it is therefore identically zero. 

Finally, it u is the weak solution of (l92l) and w = J^u, then w G C^([0, oo); L^(fi)) 
and it satisfies ^ 

Aw{t) = w{t) - [ /(r)dr Vt. 







(This is an equality as elements of H ^(f2) for all t.) However, the right hand side of 
the latter expression is a continuous function with values in L^(f2) and therefore w G 
C([0,oo);D(A)). □ 

A. 5 A simple generalization 

Consider now a closed subspace V such that Hq{Q) C C H^{Q) and that V does not 
contain non-zero constant functions, so that there exists Go such that \\u\\q < Co||Vu|| 
for all u &V. We then consider the set 

D := {ueV : Aue L^(fi), (Vm, Vf )q + (Am, f = Mv e H\n)}, 

endowed with the norm ||A ■ We can thus obtain a complete orthonormal set of 
eigenf unctions 

(pn e D, -A(j)n = Xn4>n- 

The entire theory can be repeated for these more general boundary conditions, substitut- 
ing Hq{Q) by V, -D(A) by D and H~^{Q) by the representation of V that arises from 
identifying L'^{Q) with its dual. In this case A : V ^ V is not the distributional Lapla- 
cian since elements of V cannot be understood as distributions unless V = Hq{Q). In 
any case, the results of Propositions IA.5| IA.6I and I A. 71 can be easily adapted to this new 
situation, namely. Proposition 18.31 is just a particular case. 
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